function C = convolve_monte(G,z,r,R,S0)
C=0;
[m,n]=size(G);
locz=ceil(z/0.01);
if(locz==0)
    locz=1;
end
Green=G(locz,:);
% Iphi=thetaint(r,rpp,R);
for(i=1:1:n)
    rpp=((i-1)*0.01+i*0.01)/2;
    C=C+2*pi*S0*G(locz,i)*thetaint(r,rpp,R)*rpp*0.01;
end
    






